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ABSTRACT 


The  influence  of  density  fluctuations  on  both  interface  and  volume 
reverberation  will  be  numerically  examined  in  this  work.  Using  the  same  reverberation 
geometry  and  environmental  parameters  as  defined  in  previous  works, several 
numerical  analyses  will  be  conducted  for  continuous  wave  (CW)  signal  to  predict  mean 
reverberation  structures  and  for  broadband  pulse  signals  to  generate  complex 
reverberation  structures  in  the  time-domain.  The  reverberation  model  is  based  on  the 
parabolic  equation  (PE)  approximation.  Scattering  is  assumed  to  be  dominated  by  small- 
scale  Bragg  scatter  while  the  propagation  modeling,  based  on  a  well-documented  PE 
model,  incorporates  multipath  effects  due  to  larger  range-dependent  structures.  The 
incorporation  of  density  fluctuations  in  the  PE  model  is  a  new  approach.  It  was  observed 
that  the  influence  of  bottom  volume  density  perturbations  is  to  reduce  later  (long-range) 
levels  relative  to  earlier  levels  but  does  not  appreciably  affect  the  structure.  It  was  also 
noted  that  the  CW  analysis  is  unable  to  capture  coherent  structure  of  volume  RPE  due  to 
inability  to  resolve  multipath  influence.  Therefore,  the  vertical  correlation  analysis  is  only 
valid  for  broadband  pulse  calculations.  Erom  broadband  correlation  calculations,  the 
volume  reverberation  may  decorrelate  more  rapidly  than  interface  reverberation. 
Additionally,  spectral  analysis  of  the  signals  suggested  that  response  of  interface 
reverberation  is  flatter  with  a  slope  on  the  order  of -0.125  for  both  CW  and  broadband 
data.  Volume  reverberation  displayed  a  -0.75  slope  for  CW  and  -0.25  slope  for 
broadband. 
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I.  INTRODUCTION 


The  influence  of  the  scattering  of  acoustic  energy  from  the  ocean  bottom  and  sub¬ 
bottom  has  always  been  important  to  the  Navy  because  of  reverberation  limitations  in 
active  sonar  systems.  The  main  causes  of  the  scattering  of  acoustic  energy  are  the 
impedance  contrast  at  the  ocean  bottom,  associated  rough  interface  scales,  and 
inhomogeneities  within  the  sediment  volume  producing  impedance  contrasts  of  a  separate 
spatial  structure.. 

In  the  early  90 ’s,  the  Office  of  Naval  Research  (ONR)  sponsored  a  special 
research  program  to  study  deep  ocean  reverberation  known  as  the  Acoustic  Reverberation 
Special  Research  Project  (ARSRP).'-^^  Based  on  the  results  obtained  from  this  set  of 
experiments,  the  dominant  feature  contributing  to  high  reverberation  level  was  the 
interaction  between  the  propagation  and  the  bottom  interface  topography.  ONR  has 
recently  sponsored  another  reverberation  study  within  the  ASIAEX  program  which  will 
record  reverberation  signals  in  the  shallow  water  littoral  zones  of  the  East  China  Sea.  In 
comparison  to  the  previous  study  of  the  deep  ocean,  this  area  contains  a  much  smoother 
bottom  interface  as  well  as  a  more  penetrable  sediment  interface. 

To  a  large  degree,  the  general  structure  of  the  reverberation  coincided  well  with 
the  predicted  two-way  transmission  loss.'-^^  Hence,  the  signals  measured  may  be 
predictable  via  propagation  modeling.  In  the  previous  work,  the  propagation  model, 
known  as  the  Monterey-Miami  Parabolic  Equation  Model  (MMPE),  has  been  developed 
to  incorporate  the  effects  of  bottom  interface  roughness  and  volumetric  sound  speed 
perturbations.  In  addition,  the  post-processing  routines  in  MATLAB  were  established  to 
compute  the  reverberation  loss  and  other  signal  analyses  on  the  computed  reverberation. 
The  spectral  nature  of  the  reverberation  signals  was  also  investigated. 

When  comparing  the  perturbed  data  with  unperturbed  data  in  the  environment, 
several  analyses  were  performed  for  both  CW  and  broadband  signals.  The  results  were 
summarized  as  follows:'^'^ 

•  Reverberation  level  consistently  showed  “201og  r”  drop-offs  (two-way 
cylindrical  spreading)  regardless  of  the  presence  of  perturbation. 
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•  Correlation  across  depth  showed  that  the  interface  reverberation  correlated 
better  than  the  volume.  This  also  applied  to  the  case  without  perturbations. 
Thus,  the  perturbations  were  not  significant  factors  to  affect  correlation. 

•  In  the  analysis  of  power  spectral  density  (PSD),  the  volume  reverberation 
showed  less  drop-off  for  the  perturbed  data.  However,  both  interface  and 
volume-perturbed  plots  showed  small  and  large-scale  deviations  from  the 
unperturbed  plots.  Broadband  analysis  showed  that  repeating  structures  were 
preserved  while  shorter  scale  perturbations  may  have  augmented  the  spectra. 

•  For  the  power  ratio  spectral  density  (PRSD)  analysis,  it  displayed  a  rather  flat 
spectrum  for  the  interface  reverberation  data  while  volume  reverberation 
showed  roll-off  at  higher  wavenumbers.  This  appeared  to  be  consistent  in  both 
the  CW  and  broadband  analyses. 

•  Through  signal  processing,  the  volume  perturbation  spectrum  appeared  more 
amenable  to  extraction.  However,  the  impact  of  the  interface  perturbations 
could  not  be  concluded  at  that  time  since  the  spectral  analyses  did  not  provide 
any  observable  characteristics. 

In  considering  the  effects  of  scattering  from  the  ocean  bottom,  impacts  of  the 
bottom  interface  roughness  and  sediment  sound  speed  perturbations  have  been 
incorporated  into  the  previous  model.  Based  on  the  same  reverberation  geometry  and 
environmental  parameters  as  defined  in  the  past  work,  the  objective  of  the  present  work  is 
to  investigate  the  influence  of  density  fluctuations  by  utilizing  the  improved  propagation 
model  named  as  MMPEREVDENS,  and  to  analyze  the  output  data  from  the  model  by 
implementing  the  post-processing  and  signal  processing  techniques  in  MATLAB. 
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II.  NUMERICAL  METHODS  AND  IMPLEMENTATION 


A,  REVERBERATION  THEORY 

The  sea  eontains  an  expansive  volume  of  aeoustie  scatterers  of  all  sizes  and 
interaetions,  from  interfaee  to  volume  seattering,  from  small  features  produeing  Bragg 
seatter  to  large  features  such  as  sea  mounts  and  pinnacles  producing  mostly  specular 
reflection  of  larger  wavelengths.  The  sum  total  of  the  scattering  contributions  from  all 
scatters  is  called  reverberationP^  The  definition  of  acoustic  reverberation  is  arguably  the 
received  acoustic  energy  that  hasn’t  reflected  specularly  from  the  interfaces  or  been 
refracted  by  volume  features.  In  this  thesis,  two  types  of  reverberation  are  studied, 
water/bottom  interface  scattering  and  sub-interface  volume  fluctuation  scattering.  The 
primary  distinguishing  features  between  interface  and  volume  reverberation  are  the 
impedance  contrasts  across  the  boundaries  and  the  wavelength  or  pulse  length  scales 
relative  to  the  scattering  volume.  Since  much  of  the  scattering  mechanisms  work  the 
same  whether  bistatic  or  monostatic,  only  monostatic  reverberation  will  be  studied  since 
its  numerical  implementation  is  much  easier. 

In  decibel  units  relative  to  1  |iPa  and  reference  length  scale  of  1  m,  we  define  the 
mean  reverberation  pressure  level,  RPL,  as 

i?PT  =  101og  (W\  =SL  +  DI^+DI^+mog  ^  -Rhv  (2-1) 

\Prefl\  L^oJ 

where  SL  is  the  source  level,  are  the  directivity  indices  for  the  transmitter/receiver, 

A/4^  is  the  ensonified  area  in  the  horizontal  direction,  R^^  is  the  reference  distance,  and 
RL^  ^  are  defined  as  the  reverberation  loss  for  either  the  bottom  or  the  volume. 

1.  Bottom  Interface  Scattering 

The  reverberation  loss  for  the  bottom  interface  is  defined  by^^^ 

RL,=2TL,-S,  (2.2) 

and 
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5,=101og 


(2.3) 


V  <^0  y 


K{2k,) 


where  7Z^  is  the  average  transmission  loss  (based  on  long  wavelength  components)  from 
source  to  the  scattering  patch  at  the  bottom,  is  the  full-wave  scattering  strength  due  to 
the  small-scale  interface  roughness,  is  the  wavenumber,  Cq  is  a  reference  sound  speed 
and  JV^ilkg)  is  the  2-D  spectrum  of  the  interface  roughness  evaluated  at  the  Bragg 

wavenumber. 

2,  Volume  Scattering 

For  the  volume  scattering,  the  volume  reverberation  loss  cannot  be  expressed 
simply  in  terms  of  the  two-way  transmission  loss  but  instead  must  be  the  integral  over 

depth  of  the  quantity  at  each  range  r,  where  n{r,z)  is  the  approximate 

refractive  index  based  on  only  long  wavelength  perturbations,  and  Y{r,z^  is  the  field 
function  of  the  two-way  propagation.  The  reverberation  loss  for  the  volume  is  then 
defined  by'-^^ 


i?I,=-201og 


1 

r 


J  h{r,z)\y/{r,z)'\dz 


(2.4) 


and 


‘S,=101og 


(2.5) 


In  the  above  expression,  is  the  volume  scattering  strength,  JV2^(2kg)  is  the  2-D 
horizontal  spectrum  of  the  volume  fluctuations,  which  is  assumed  isotropic  and 
independent  of  depth  evaluated  at  the  Bragg  wavenumber. 
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B,  MONTEREY-MIAMI  PARABOLIC  EQUATION  (MMPE)  MODEL 

The  parabolic  equation  (PE)  method  was  introduced  into  underwater  acoustics  in 
the  early  1970’s  by  Tappert.'-"^^  It  is  a  numerical  model  to  solve  acoustic  wave  propagation 
problems  in  the  ocean.  The  PE  method  has  become  the  most  popular  wave-theory 
technique  for  solving  the  underwater  acoustic  range-dependent  problem.'-^^  The  MMPE 
Model  is  based  upon  the  parabolic  approximation  and  thus  a  brief  description  of  this 
approach  is  useful. 

If  the  monofrequency  acoustic  pressure  field  is  written  in  cylindrical  coordinates 
r  =  {r,(p,z)  as 

P{r,(p,z,0)t)  =  p{r,(p,z)e~‘“  (2.6) 


and  substituted  into  the  wave  equation,  we  obtain  the  Helmholtz  equation 

VV(r)  +  ^;7(r)  =  0  (2.7) 

c 

where 


1  a  a  1  a'  a^ 

- v - 1 - 1 - , 

r  dr  dr  r^  d(p  dz^ 


(2.8) 


The  ocean  environment,  and  thus  the  underwater  acoustic  field,  tends  to  exhibit 
only  weak  azimuthal  dependence.  This  allows  us  to  ignore  the  azimuthal  terms,  and 
define  the  complex  acoustic  pressure'-^^ 

P[r,z,cot)  =  p[r,z)e~““ .  (2.9) 


To  simplify  the  Helmholtz  equation  and  account  for  the  cylindrical  spreading,  we 
define'^^^ 

p{r,z)  =  ^u{r,z)  (2.10) 

The  1  /  4r  term  accounts  for  the  azimuthal  spreading  and  u{r,z)  is  the  two-dimensional 
pressure  field.  Substituting  Eq.  (2.10)  into  Eq.  (2.7)  yields^^^ 
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-+- 


dr  dz^ 


■+kl 


+ 


M  =  0  . 


(2.11) 


This  is  referred  to  as  the  uneoupled  azimuth  approximation.  The  Helmholtz  equation  ean 
be  faetored  by  introdueing  the  operator  notation 


P  =  — 
dr 


(2.12) 


and 


Q^^  =  iju+£+iy 


(2.13) 


where 


£  =  n^  -  \  and  jU  = 


1  a' 

dz^  ■ 


(2.14) 


Rewriting  Eq.  (2.1 1),  we  obtain 

(2-15) 

The  proper  faetorization  of  the  outward  propagating  field  is  obtained  by  defining'-^^ 

u  =  ar'E  (2.16) 

then  the  outgoing  wave  must  satisfy'-^^ 

;)vi/ 

(2.17) 


All  PE  models  are  based  on  Eq.  (2.17).  It  represents  the  complete  description  of  the 
forward  propagating  acoustic  energy  in  the  waveguide.  The  difference  among  each  model 
is  the  approximations  of  the  operator  Qop  and  the  method  of  generating  solutions  to  this 
equation. 

To  develop  a  numeric  algorithm  for  solving  the  parabolic  equation,  the  acoustic 
field  may  be  decomposed  into  a  slowly  modulating  envelope  function  and  a  rapidly 
varying  phase  term  that  oscillates  at  the  acoustic  frequency.  The  envelope  function,  or 
field  function  y/(r,z),  is  defined'^*’^ 
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(2.18) 


or,  in  terms  of  the  aeoustie  pressure, 

.  (2.19) 

V  r 

The  parabolie  equation  for  the  field  funetion  is  then  defined 

^  =  -KW + iKQcpW  =  (2.20) 

where 

Kp='^-Qop  (2.21) 

is  a  Hamiltonian-like  operator  whieh  defines  the  evolution  of  the  PE  field  funetion  in 
range. 

The  relationship  between  values  of  ^^at  different  ranges  ean  be  expressed  as'-^^ 

^(r-l-Ar)  =  d>(r)^(r)  (2.22) 


where  ^>(r)  is  a  propagator  that  marehes  the  solution  out  in  range.  The  MMPE  model 
employs  a  split-step  Eourier  (PE/SSE)  method'-^^  to  provide  a  representation  of  the 
propagator  ^>(r).  This  method  is  utilized  primarily  based  on  the  speed  and  simplieity  of 
the  PE/SSE. 

It  is  essential  to  the  SSE  algorithm  that  the  different  terms  within  Hop  operator 
must  be  separated  whieh  requires  an  approximation  to  the  square  root  operator,  .  The 

approximation  of  the  Hamiltonian  operator  used  in  the  MMPE  eorresponds  to  the  so- 
ealled  wide-angle  approximation  (WAPE)'^^^ 

H  +U  (2.23) 

where 


Top{k)  =  \- 


1 

2kl  dz^ 


(2.24) 
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and 


t',,  =-['<-!]  ■ 


(2.25) 


The  propagator  function  4>(r)  can  then  be  expressed  as'^' 


d>(r) = 


Q  I  ^  ^  ^  ^ 


(2.26) 


The  PE/SSF  algorithm  employs  the  operator  e  in  the  ^z-domain  to  simplify  the 
calculations  where'-^^ 


T  =1-  1 - ^ 

u. 


(2.27) 


With  the  fast  Fourier  transform  employing  the  convention 


yr{z)  =  FFT[^{K)) 


(2.28) 


Y[K)  =  IFFT(yf[z)) 


(2.29) 


the  PF/SSF  implementation  can  be  represented  by'^' 


\y  I  !  \y 

/  \  — ^^(r+Ar.z)  -it  Art  (t'i  ~‘K — V^Ar,z) 

y/[r  +  Ar,z)  =  e  ^  "  xFFT \e  "^^^^IFFT  e  ^  "  : 


xy/{r,z)  .  (2.30) 
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C.  IMPLEMENTATION  OF  THE  REVERBERATION  PROBLEM  IN  MMPE 

Having  briefly  described  concepts  of  reverberation  theory  and  the  MMPE  model, 
in  this  section  we  will  concentrate  on  the  theoretical  descriptions  for  generating 
perturbations  to  both  interface  and  volume,  including  the  influence  of  volume  density 
fluctuations.  The  incorporation  of  these  effects  into  the  MMPE  model  is  also  discussed. 
While  the  theoretical  basis  for  modeling  the  interface  roughness  is  based  on  the  work  of 
Goff  and  Jordon, the  development  for  the  volume  is  based  on  Yamamoto’s  work.'-^*’^ 

1.  Interface  Roughness 

We  assume  a  2-D  (two  dimensional)  interface  spectrum  of  the  form'-^^ 

W^{k^)  = - ^ - j-  and  =  ^Jk^  +  L^  (2.31) 


where  k^  is  the  horizontal  spatial  wavenumber  vector,  ^  is  a  normalization  factor, 
is  a  correlation  length  scale,  ^  is  the  spectral  exponent  and  K  and  L  are  the  horizontal 
wavenumbers  in  the  x-  and  y-directions,  respectively. 

If  the  2-D  spectrum  {k^)  is  assumed  to  be  independent  of  direction  (isotropic), 
then  the  normalization  factor  //  can  be  defined  in  terms  of  the  root-mean-square  (rms) 
roughness  (J^  by  requiring'-^^ 


2;rJ  W^{k^)k^dk^  = 

0 


(2.32) 


which  leads  to 


/^  = 


U 


K 


<7  L 


2 

'corr  ' 


(2.33) 


We  may  simply  evaluate  W^ik^)  at  k^  =  Ik^  for  the  scattering  amplitude  since  it 
is  caused  by  Bragg  scatter  (evaluated  along  the  line  of  propagation  for  monostatic 
reverberation).  However,  for  the  long-wavelength  interface  roughness,  we  need  to 
evaluate  the  full  spectrum  effect.  This  could  be  done  by  determining  the  1-D  (one 
dimensional)  spectrum  along  a  slice  and  simply  use  a  1-D  transform  which  is  defined  as'-^^ 
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W,{K)=\  W^{K,L)dL 


(2.34) 


Or  in  cylindrical  coordinates,  we  have'-^^ 


=  4  /  ,  'V,sk,)dk,  = 


where 


-1 


r=- 


2  2 


A 


J 


(2.35) 


(2.36) 


In  order  to  generate  a  1-D  roughness  realization  from  Eq.  (2.34)  or  Eq.  (2.35),  we 
transform  the  1-D  amplitude  speetrum  that  has  been  sealed  by  a  random  amplitude  and 
phase.  In  other  words,  we  define  the  realization  as'-^^ 


7(x)=  j  S,{K)e'^dK 


(2.37) 


where 

S,{K)  =  [W,{K)Y'^  (2.38) 

and  A  and  ^are  random  numbers  for  all  values  of  K. 

The  random  phase  and  the  amplitude  of  eaeh  eomponent  can  be  obtained  from  Eq. 
(2.39)  and  Eq.  (2.40),  respeetively,'-'^ 

0  =  27rr^  (2.39) 

A  =  ^-ln{r,),  (2.40) 

where  both  ri  and  r2  are  independent  uniformly  distributed  random  variables  in  the 
interval  [0,1].  In  praetiee,  we  eould  simply  use'^'^ 

If](^)  =  (l  +  EL^^p"7.  (2.41) 
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2,  Volume  Sound  Speed  Fluctuations 

The  sediment  volume  sound  speed  perturbations  may  be  modeled  by  a  3-D  (three 
dimensional)  volume  spectrum  given  by'-^^ 

W^{K,L,M)  =  A" (k^  +  L^)  +  M^y  ^ (2.42) 

where  B  is  the  spectral  strength  constant,  P  is  the  spectral  exponent,  and  A  =  —  =  —  is 

d-^  Cl'2^ 

the  horizontal-to-vertical  aspect  ratio  describing  the  anisotropy  of  fluctuations  in  the 
sediment. 

To  evaluate  the  reverberation  due  to  the  volume  perturbations,  we  need  an 
expression  for  the  2-D  horizontal  spectrum  (assuming  strongest  scattering  near 
horizontal).  It  is  defined  as'-'^ 

W^{K,L)  =  j  W^{K,L,M)dM  .  (2.43) 

Substituting  Eq.(2.42)  into  Eq.(2.43),  we  obtain 

2  P 

=  +  L^)  +  M^y^~'  dM  .  (2.44) 

^  0 

Eor  /3  =  2  ,  Eq.  (2.44)  can  be  reduced  to 

W^{K,L)  =  ^^A^(k^ .  (2.45) 

Eor  the  values  of  5  ~  5x10“"^ and  A  ~  5  chosen  from  Yamamoto’s  findings'-^°\  Eq.  (2.45) 
reduces  to'-'^ 

_  3 

W^{K,L)  =  a[K^  +  L^)~  ^  (2.46) 

where 
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(2.47) 


a  =  — =  5x10"' 

2A 

However,  to  simplify  the  forward  propagation,  we  need  only  the  2-D  vertical 
speetrum  in  the  (r,z)  plane.  It  is  defined  by'-^^ 

W^{K,M)=  \w,{K,L,M)dL  =  ^ - j  +  f')  +  M']  ^  (2.48) 

For  P-2,  Eq.  (2.48)  becomes 

_  3 

W^{K,M)  =  a\25K^ (2.49) 

where 

cr'  =  1.25x10"'  .  (2.50) 

To  generate  2-D  vertieal  volume  sound  speed  fluetuation  realizations,  a 

realization  ean  be  defined  as'-'^ 

c,S{x,z)  =  \\s^{K,M)e'^e^^dKdM  (2.51) 

where 

S^{K,M)  =  [W^{K,M)Y^  ■  (2.52) 

Notiee  that  we  really  are  generating  a  series  of  vertical  realizations  at  eaeh  range  step, 

sinee  we  have  treated  sound  speed  perturbations  in  the  volume  in  the  vertical. 

Similar  to  the  interfaee,  the  2-D  random  phase  and  the  amplitude  variations  ean  be 
obtained  in  Eq.  (2.53)  and  Eq.  (2.54),  respeetively,^'^ 

9{K,M)  =  27rr^{K,M)  (2.53) 

A{K,M)  =  ^-{n{r,{K,M))  (2.54) 

where  both  ry{K,M)  and  r2{K,M)  are  now  a  matrix  of  uniformly  distributed  random 


numbers  in  [0,1]. 

In  praetiee,  we  use^^^ 
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^  - 1 
2 


(2.55) 


Eqs.  (2.41)  and  (2.55)  are  the  results  for  the  interface  and  volume  perturbations, 
respectively.  These  are  the  generic  spectral  models  used  in  generating  the  realizations 
first  in  MATLAB  and  then  for  implementation  in  the  MMPE  model. 

3,  Density  Fluctuations  in  Sediment 

Variability  in  density,  p,  is  incorporated  into  the  PE  model  by  defining  the 
effective  index  of  refraction'-^ 


/2 


2^  1 

n  H - ; 

2^ 


1  V72  3 

—  VV  — 
P  2 


1 


,  V/7 

VP  J 


(2.56) 


Consistent  with  the  numerical  treatment  that  assumes  the  environment  is  range- 
independent  over  a  range  step,  and  the  fact  that  sediment  properties  are  largely 
horizontally  stratified,  we  may  simplify  this  to''^" 


'2 


■  n  -I-- 


2kl 


1  av 

p  dz^ 


yP 


(2.57) 


Eor  the  forward  problem,  the  sound  speed  index  of  refraction  is  based  only  on  large  scale 
features,  such  that'^'-' 


'2  y-  \  ‘-0 


(2.58) 


and 


(2.59) 


where  is  the  mean  bottom  sound  speed  at  the  interface,  =  y  is  normalized 

gradient  of  bottom  sound  speed,  and  <5)  is  the  zero-mean  random  perturbation  for  the 
long  wavelength  component. 

According  to  the  analysis  of  Yamamoto  (1996),  the  relative  fluctuations  in  density 
relate  to  the  relative  fluctuations  in  sound  speed  according  to'^'^ 
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(2.60) 


^  ^  2(P, -Po)  Sc  ^  2^^ 
Po  ~  Pr  ^0  ^0 


where 


z,„(l  +  ^z)  and  r-- 


Pr-P, 

^Po-Pr 


(2.61) 


Notiee  that  and  Cq  are  the  averaged  values  in  the  sediment  and  =  2650^y  3  is  the 
density  of  the  grain.  We  may  then  write'-^'^ 


p  =  p^  +  Sp  =  p,  \  +  -^dc  . 


(2.62) 


Taking  the  first  and  seeond  partial  derivatives  of  Eq.  (2.63)  with  respeet  to  depth,  z,  and 
negleeting  depth  gradients  in  either  CgOr  p^ ,  we  obtain'^' 


dz  Ca  dz 


(2.63) 


(2.64) 


With  the  sound  speed  fluetuations  being  defined  by^ 


=  CgS  (x,  z)  =  {K,M)e‘^e^"dKdM 


(2.65) 


where 


^2  =  [1F2'(.^,M)]^ 


W(K,M) 


(2.66) 


Substituting  Eq.  (2.65)  into  Eq.  (2.63)  and  Eq.  (2.64)  yields'^ 


^  JJ {iM)S^  {K,M)e‘^e‘^^dKdM 


(2.67) 
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(2.68) 


az  Cg 

The  sediment  effeetive  index  of  refraetion  ean  now  be  derived  by  substituting  Eq.  (2.67) 
and  Eq.  (2.68)  into  Eq.  (2.57)  whieh  beeomes^'^^ 

nl^{x,z)  =  nl{x,z)  +  ^\_a{x,z)  +  /3{x,z)~]  (2.69) 

where 

a{x,z)=  \\m^S^  {K,M)e'‘^e^^dKdM  (2.70) 

c,[x,z)^^ 

and 

r 

/5[x,z)  =  +  y.  l\MS^[K,M)e'’^d^^dKdM  .  (2.71) 

These  parameters  will  be  eomputed  in  parallel  with  in  MMPEREV.  In  the  sediment, 
there  is  then  the  additional  propagator  term'-^'^ 

=  (2.72) 

where 

=  +  (2-73) 

HKq 
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D,  TIME-DOMAIN  PROCESSING 

Previously  discussion  of  the  reverberation  loss  for  either  the  bottom  or  the 
volume,  ,  is  based  on  continuous  wave  (CW)  analysis.  In  order  to  obtain  the 

structure  of  the  pulse  propagation  in  time,  and  more  properly  treat  coherent  interference 
effects  by  separating  multipaths,  we  need  to  implement  the  MMPE  model  over  a 
spectrum  of  frequencies.  In  the  following  discussion,  the  time-domain  analysis  of  the 
interface  and  the  volume  will  give  us  the  general  picture  of  the  two-way  travel  time 
structure  of  the  reverberation  loss.  We  can  then  arrive  at  the  reverberant  field  at  each 
range  step  and  continue  the  propagation  through  the  entire  water  column  of  interest. 

1.  Time-Domain  Analysis  of  the  Interface 

The  geometry  of  a  two-way  return  from  a  scattering  patch  adapted  from  Smith 
and  Cushman'-^^  is  shown  in  Figure  1. 


The  travel  time  of  the  two-way  pressure  field  at  the  receiver  is  the  convolution  of 
two,  one-way  fields  in  the  time-domain'-^^ 

Pl-«ay,b  .  'Z')  =  J  Pyn  ('h  .  t)PyRh  (fn  ^  'Z'  “  t)dt  (2.74) 

where  and  are  the  forward  propagated  pressure  fields  from  the  transmitter  and 

receiver  to  the  scattering  point  evaluated  at  the  bottom  interface,  respectively.  Notice  that 
the  receiver  and  transmitter  are  not  necessary  to  be  co-located  in  the  water  column  and  by 
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reciprocity,  the  propagated  field  from  the  receiver  to  the  scattering  point,  is  the 


same  as  the  propagated  field  from  the  scattering  point  to  the  receiver 


Because  the  time-domain  convolution  of  the  two  field  functions  is  also  the  scalar 
multiplication  of  these  functions  in  the  frequency  domain,  the  two-way  field  in  the 
frequency  domain  from  the  interface  can  be  expressed  as'-^^ 

P2-.ay,b(^mJ)  =  P  f)  P  yRhir^J)  (2.75) 

where 


Pynir^J) 


(2.76) 


and 


PyRbir^J) 


(2.77) 


Furthermore,  the  two-way  travel  time  structure  of  the  reverberation  loss  for  the  bottom 
interface,  RL^ ,  due  to  a  single  bottom  patch  can  be  derived  as'-^^ 

Pl-y.ay,b  (f„  ’  0  =  ^  J  Pl-y,ay,b  (fn  ^  df  (2-78) 

where  ^4  is  a  constant  which  contains  all  the  other  factors  needed  to  define  reverberation 
loss. 

2,  Time-Domain  Analysis  of  the  Volume 

By  applying  a  Fourier  transform  to  the  time-domain,  the  reverberant  field  due  to 
each  depth/range  point  is  defined  by'^'^ 

P2-^ay,v  =  «(f«  ^  ^)Pyr  (f«  ^  f)PyR  f)  (2-79) 

where  now  the  two-way  reverberation  signal  is  computed  for  every  grid  point  of  interest 
( z  >  always)  at  a  particular  frequency,  f,  and  n{r^ ,  z)  is  the  local  index  of  refraction  at 
the  grid  point.  In  order  to  provide  the  same  weighting  used  in  the  CW  analysis,  the 
multiplication  term,  n{r^,z) ,  is  required  in  Eq.  (2.79). 
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Additionally,  the  two-way  travel  time  structure  of  the  reverberation  loss  for  the 
volume,  RLv ,  can  be  derived  as'-'^ 

P2-^ay,v  =  ,  t)dz  (2.80) 

Z>Zj 

Note  that  the  constant  A  and  B  in  Eq.  (2.78)  and  Eq.  (2.80)  contain  all  the  other  factors 
needed  to  define  reverberation  loss. 

The  implementation  of  the  above  discussions  for  both  the  interface  and  volume 
will  be  discussed  in  the  post-processing  analysis  in  Chapter  III. 
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III.  RESULTS 


A,  DESCRIPTION  OF  EXPERIMENTS 
1.  Reverberation  Geometry 

In  experiments,  an  vertieal  line  array  (VLA)  with  16-element  was  ehosen  to 
provide  the  reverberation  measurements/eomputations  for  both  monostatie  (eoineident) 
and  bistatie  (vertieal  separated)  eases.  The  array  was  loeated  vertieally  in  the  water 
eolumn  and  eaeh  element  of  the  VLA  was  assumed  to  be  a  point  souree.  The  geometry  of 
VLA  and  seattering  pateh  was  shown  in  Figure  2. 

From  Figure  2,  we  observed  that  the  water  eolumn  was  spanned  by  16-element 
vertieal  line  array  from  20  m  to  80  m  with  4  m  apart  in  depth  between  eaeh  element.  The 
souree  was  loeated  at  depth  of  48  m  with  all  16  elements  reeeiving  the  reverberation.  The 
range  r  is  the  mean  horizontal  distanee  of  the  VLA  to  a  small  seattering  pateh. 


Figure  2.  Geometry  of  VLA  and  Seattering  Pateh  [From  Ref.  1] 
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2.  The  Environmental  Parameters 

The  environmental  parameters  used  in  the  MMPE  model  were  speeified  below: 


Filename/Parameter 

Value 

Remarks 

Main  Control  File:  pefiles.inp 

Number  of  depth  points 

256 

Radix-2  integer  required 

Minimum  depth 

Om 

Maximum  depth 

400  m 

Number  of  range  steps 

833 

Minimum  range 

Om 

Maximum  range 

5.0  km 

Range  step  size 

6  m 

Maximum  computed  depth 

400  m 

Reference  sound  speed 

1500  m/s 

Source  File:  pesrc.inp 

Source  depths 

Varying 

Array  elements  at  20,  24,  28,  32,  36, 

40,  44,  48,  52,  56,  60,  64,  68,  72, 

76  and  80  m  depths. 

Center  frequency 

250  Hz 

Frequency  bandwidth 

250  Hz 

No.  of  Frequencies 

512 

Radix-2  integer  required 

Sound  Speed  File:  pessp.inp 

Water  column  sound  speed 

1500  m/s 

Range  independent 

No.  ofSSPs 

1 

Bathymetry:  pebath.inp 

Mean  bottom  depth 

100  m 

Range  independent 

No.  of  depth  points 

1 

Bottom  properties:  pebotprop.inp 

Bottom  sound  speed 

1700  m/s 

Sound  speed  gradient 

0 

Relative  density 

1.0 

No  density  variations 

Compressional  attenuation 

0.2  dB/km/Hz 

Shear  speed 

0 

Not  modeled 

Shear  attenuation 

0 

Not  modeled 

Deep  Bottom  Bathymetry:  pedbath.inp 

Depth 

3000  m 

Deep  Bottom  Properties:  pedbotprop.inp 

Deep  bottom  sound  speed 

2000  m/s 

Sound  speed  gradient 

0 

Relative  density 

3.0 

No  density  variations 

Compressional  attenuation 

0.25  dB/km/Hz 

Shear  speed 

0 

Not  modeled 

Shear  attenuation 

0 

Not  modeled 

RMS  Perturbations  (input  to  MMPEREV  during  run) 

Interface  roughness 

1  m 

Volume  sound  speed  fluctuation 

15  m/s 

Table  1.  MMPEREV  Input  Environmental  Properties  [After  Ref  1] 
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Few  points  should  be  noted: 

•  Seven  control  and  environmental  input  files  were  required  by  MMPE  to 
perform  the  computations. 

•  Most  of  the  parameters  used  for  the  bottom  and  deep-bottom  properties  were 
assumed. 

•  To  perform  broadband  as  well  as  CW  analysis,  a  center  frequency  of  250  Hz 
was  chosen  with  a  250  Hz  bandwidth  divided  into  512  discrete  propagation 
frequencies. 

•  In  a  shallow  water  environment,  the  mean  bottom  depth  was  assumed  to  be 
100  m  with  no  mean  slope. 

•  Both  the  number  of  frequencies  and  depth  points  have  to  be  radix-2  integer 
because  of  the  Fast-Fourier  transform  (FFT)  that  was  used  to  compute  spectral 
components  in  MMPE  model. 
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B,  POST-PROCESSING  1-  REVERBERATION  LOSS  (RL) 

IN  CONTINUOUS  WAVE  (CW)  AND  TIME-DOMAIN  (BB) 

1,  CW  Reverberation  Analysis 

Equations  needed  to  formulate  the  MATLAB  implementation  for  the  interface 
and  the  volume  reverberation  loss  were  stated  as  follows: 

•  For  the  interface  reverberation  loss, 


RLf^  (r)  =  y4-201og 


r 


(3.1) 


where  A  is  a  constant  which  accounts  for  the  other  parameters  in  Eq.  (2.2). 
•  For  the  volume  reverberation  loss,  RL^ 


i?4(r)  =  5-201og 


h[r,z)  y/j.[r,z)  l//^^[r,z)dz 


(3.2) 


where  B  is  another  constant  which  accounts  for  the  other  terms  in  Eq.  (2.4). 

2,  CW  Analysis  Results 

With  the  constants  A  in  Eq.  (3.1)  and  B  in  Eq.  (3.2)  set  equal  to  zero,  a  source 
depth  of  48  m,  receiver  depth  of  40  m  and  frequency  of  250  Hz,  the  interface  and  volume 
reverberation  loss  for  perturbed  data  with  and  without  density  fluctuations  were  shown  in 
Figures  3  and  4,  respectively.  We  observed  that  the  reverberation  structure  of  interface 
and  volume  for  perturbed  data  with  density  fluctuations  displayed  a  higher  reverberation 
loss  than  perturbed  data  without  density  fluctuations.  The  differences  of  these  two  were 
less  than  10  dB. 

Eikewise,  the  colormaps  of  receiver  depths  versus  ranges  for  both  the  interface 
and  volume  reverberation  were  used  to  analyze  the  effects  of  density  fluctuations.  From 
Figures  5  and  6  (interface  reverberation)  and  Figures  7  and  8  (volume  reverberation),  we 
can  conclude  that  perturbed  data  with  and  without  density  fluctuations  showed  very  small 
differences. 
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Source  depth:48m,  Receiver  depth:40m,  Frequency:250Hz 


Figure  3.  Interface  Reverberation  Loss  for  Two-Way  Transmission  Vs  Range 


Figure  4.  Volume  Reverberation  Loss  for  Two-Way  Transmission  Vs  Range 
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Figure  5. 


Figure  6. 


Interface  Reverberation  Loss  -  Receiver  Depth  Vs  Range  With  Perturbation  and 

With  Density  Fluctuations 


Source  Depth:  48m,  Frequency:  250Hz 
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Interface  Reverberation  Loss  -  Receiver  Depth  Vs  Range  With  Perturbation  and 
With  No  Density  Fluctuations 
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Source  Depth:  48m,  Frequency:  250Hz 


I 

Range  (km) 


Figure  7.  Volume  Reverberation  Loss  -  Receiver  Depth  Vs  Range  With  Perturbation  and 

With  Density  Fluctuations 


Source  Depth:  48m,  Frequency:  250Hz 


Range  (km) 

Figure  8.  Volume  Reverberation  Loss  -  Receiver  Depth  Vs  Range  With  Perturbation  and 

With  No  Density  Fluctuations 
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3.  Broadband  Reverberation  Analysis 

To  predict  the  influences  of  a  pulsed  signal,  the  analyses  of  the  pulse  propagation 
in  the  time-domain  for  both  the  reverberation  structure  of  the  interface  and  volume  were 
necessary.  Equations  required  to  formulate  the  MATLAB  implementation  for  the 
interface  and  the  volume  reverberation  loss  were  summarized  below: 

•  For  the  interface  reverberation  loss, 


P2-^ay,b  ('h  .  /)  =  PyTh  ('h  ^  f)PyRh  ('h  ^  /) 


(3.3) 


where 


Pywi^mJ) 


¥yTbir„J)e 


ik„r 


(3.4) 


and 


PyRhir„,f) 


¥yRbir„,f)e 


ik„r 


(3.5) 


> 


> 


P2-^ay,b  ,  0  =  J  P2-«ay,b  ('h  ^  df 

M 

P 2-way, b  (^«  )  ^  1  P 2-way, b  ) 


m—l 


•  For  the  volume  reverberation  loss,  RL^ 


[1] 


> 

> 

> 


P2-way,v  ¥m  ’  ^  J)  =  "('h  ^  ^)PyT  f  )PyR  ('h  ^  /) 

P  2-way, V  ('h  ^  0  =  J  P2-way,v  ^  df 

M 

P2  -way,v  ('.)=E  P2  -way^v  (^m  ’  ) 


(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 
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4.  Broadband  Analysis  Results 

With  a  receiver  depth  of  40  m  chosen  arbitrarily  and  a  source  depth  of  48  m,  the 
general  structures  of  interface  and  volume  reverberation  for  perturbed  data  with  and 
without  density  fluctuations  in  the  time-domain  were  shown  in  Figure  9  (interface 
reverberation)  and  in  Figure  10  (volume  reverberation).  For  both  the  interface  and 
volume  reverberation,  the  comparison  between  perturbed  data  with  and  without  density 
fluctuations  had  shown  similar  characteristics  within  1  second  of  travel  time  and  differed 
slightly  beyond  this  interval.  Overall,  the  differences  of  reverberation  loss  were  within  8 
decibels. 

For  the  colormaps  of  receiver  depths  versus  ranges  for  both  interface  (Figures  1 1 
and  12)  and  volume  reverberation  (Figures  13  and  14),  the  perturbed  data  without  density 
fluctuations  displayed  more  structure  while  the  perturbed  data  with  density  fluctuations 
revealed  little  structure. 
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Source  depth:48m,  Receiver  depth:40m 


Figure  9.  Time -Domain  Interface  Reverberation  Loss  for  Two-Way  Transmission 


Figure  10.  Time -Domain  Interface  Reverberation  Loss  for  Two-Way  Transmission 


28 


Source  Depth:  48m 


Time  (sec) 

Figure  1 1 .  Interface  Reverberation  Loss  Receiver  Depth  Vs  Time  With  Perturbation  and 

With  Density  Fluctuations 


Source  Depth:  48m 


Time  (sec) 


Figure  12.  Interface  Reverberation  Loss  Receiver  Depth  Vs  Time  With  Perturbation  and 

With  No  Density  Fluctuations 
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Source  Depth:  48m 


Time  (sec) 


Figure  13.  Volume  Reverberation  Loss  Receiver  Depth  Vs  Time  With  Perturbation  and  With 

Density  Fluctuations 


Source  Depth:  48m 


u  u.J  I  I  e.  e—J  o 

Time  (sec) 

Figure  14.  Volume  Reverberation  Loss  Receiver  Depth  Vs  Time  With  Perturbation  and  With 

No  Density  Fluctuations 
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C.  POST-PROCESSING  2-  VERTICAL  CORRELATION  &  PEAK 

CORRELATION  IN  CW  AND  BB 

1.  Vertical  Correlation  in  Range  (CW) 

A  source  depth  of  48  m  and  the  eenter  frequency  of  250  Hz  were  chosen  to 
perform  the  correlation  analysis.  The  perturbed  data  composed  with  and  without  the 
density  fluctuations  in  Figures  15  and  16  portrayed  the  vertieal  correlation  for  the 
interface  while  Figures  17  and  18  showed  the  case  of  the  volume  for  the  perturbed  data 
with  and  without  density  fluctuations. 

In  general,  the  strueture  of  perturbed  data  with  and  without  density 
fluctuations  was  similar  for  both  cases  of  the  interface  and  volume.  For  the  interface 
reverberation  range  correlation,  the  perturbed  data  with  density  fluctuations  revealed 
stronger  eorrelation  throughout  depth  for  the  range  lag  between  0.3  and  0.5  km.  This 
stronger  correlation  could  also  be  found  in  the  perturbed  data  with  density  fluctuations 
for  the  range  lag  between  2.7  and  3  km  in  the  analysis  of  volume  reverberation. 

2,  Peak  Vertical  Correlation  in  CW 

To  observe  the  rate  of  change  for  the  signal  decorrelated  over  depth,  the 
values  of  peak  eorrelation  were  extraeted  from  the  above  vertical  correlation.  Figure  19 
showed  a  eomparison  of  the  peak  vertical  correlations  between  perturbed  data  with  and 
without  density  fluctuations  for  both  the  interface  and  volume.  Aceording  to  the  figure, 
we  ean  summarize  that  the  addition  of  density  fluctuations  did  not  appear  to  affeet  the 
vertieal  correlations  dramatically.  Additionally,  the  interface  showed  more  deeorrelation 
over  depth  than  did  the  volume  displayed  more  correlated  over  depth. 
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Source  Depth:  48m,  Frequency:  250Hz 


-0.5  -0.4  -0.3  -0.2  -0.1  0  0.1  0.2  0.3  0.4  0.5 

Range  Lag  (km) 


Figure  15.  Vertical  Correlation  of  Interface  Reverberation  Loss  in  Range  Vs  Relative  Depth 

With  Perturbation  and  With  Density  Fluctuations 


Source  Depth:  48m,  Frequency:  250Hz 
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Range  Lag  (km) 


Figure  16.  Vertical  Correlation  of  Interface  Reverberation  Loss  in  Range  Vs  Relative  Depth 

With  Perturbation  and  With  No  Density  Fluctuations 


32 


Source  Depth:  48m.  Frequency:  250Hz 


-3-2-1  0  1  2  3 

Range  Lag  (km) 


Figure  17.  Vertical  Correlation  of  Volume  Reverberation  Loss  in  Range  Vs  Relative  Depth 

With  Perturbation  and  With  Density  Fluctuations 


Source  Depth:  48m,  Frequency:  250Hz 


Range  Lag  (km) 

Figure  18.  Vertical  Correlation  of  Volume  Reverberation  Loss  in  Range  Vs  Relative  Depth 

With  Perturbation  and  With  No  Density  Fluctuations 
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Source  Depth:  48m,  Frequency:  250Hz 


Peak  Correlation  (dB) 


Figure  19.  Peak  Vertical  Correlation  Vs  Relative  Depth 
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3.  Vertical  Correlation  in  the  Time-Domain  (BB) 

For  the  broadband  signals,  the  vertieal  correlation  for  interface  and  volume 
was  computed  in  the  time-domain.  These  results  were  listed  in  Figures  20  and  21 
(vertical  temporal  correlation  of  interface  reverberation  loss  versus  relative  depth)  and  in 
Figures  22  and  23  (vertical  temporal  correlation  of  volume  reverberation  loss  versus 
relative  depth)  for  the  perturbed  data  with  and  without  density  fluctuations. 

Comparison  of  the  perturbed  data  with  and  without  density  fluctuations 
reveals  little  difference  between  them.  This  tiny  difference  applied  to  both  the  interface 
and  volume  reverberation  temporal  correlation.  Furthermore,  the  finer  scale  structures 
were  also  noted  in  both  the  interface  and  volume  reverberation  temporal  correlation. 

4.  Peak  Vertical  Correlation  in  BB 

In  Figure  24,  the  results  on  the  perturbed  data  computed  with  and  without 
density  fluctuations  of  peak  vertical  temporal  correlation  were  extracted  from  the 
previous  analysis  of  vertical  correlation  for  the  interface  and  volume  relative  to  a  source 
depth  of  48  m.  By  comparison  of  these  data  in  Figure  24,  the  influence  of  density 
fluctuations  was  minor  in  both  the  interface  and  volume  reverberation  loss. 
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Source  Depth:  48m 


Figure  20. 


Figure  21. 


Vertical  Temporal  Correlation  of  Interface  Reverberation  Loss  Vs  Relative  Depth 
With  Perturbation  and  With  Density  Fluctuations 


Source  Depth:  48m 
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Vertical  Temporal  Correlation  of  Interface  Reverberation  Loss  Vs  Relative  Depth 
With  Perturbation  and  With  No  Density  Fluctuations 
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Source  Depth:  48m 
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Figure  22.  Vertical  Temporal  Correlation  of  Volume  Reverberation  Loss  Vs  Relative  Depth 

With  Perturbation  and  With  Density  Fluctuations 


Source  Depth:  48m 
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Figure  23.  Vertical  Temporal  Correlation  of  Volume  Reverberation  Loss  Vs  Relative  Depth 

With  Perturbation  and  With  No  Density  Fluctuations 
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Source  Depth:  48m 


Figure  24.  Peak  Vertical  Temporal  Correlation  Vs  Relative  Depth 
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D,  POST-PROCESSING  3-  SPECTRAL  CHARACTERISTICS  IN  CW  AND 

BB 

In  order  to  extract  the  spectral  components  of  the  reverberation  data  for  both  the 
CW  and  broadband  signals,  the  analysis  was  performed  by  using  the  methods  of  signal 
processing,  such  as  the  power  spectral  density  and  power  ratio  spectral  density. 

1.  The  Analysis  of  Power  Spectral  Density  (PSD) 

The  interpretation  of  the  power  density  can  be  represented  as  the  correlation  of 
the  random  process  at  a  particular  frequency.  The  spectra  gave  quantitative  descriptions 
of  the  frequency  components  of  a  random  process.  In  other  words,  the  power  spectral 
density  can  be  defined  as  a  value  of  the  spectrum  at  a  given  radian  frequency. 

For  the  CW  signal,  the  plots  of  normalized  power  versus  wave  number  for  both 
interface  and  volume  reverberation  in  cases  of  perturbed  data  with  and  without  density 
fluctuations  can  be  shown  in  Figures  25  and  26,  respectively.  Comparison  of  these 
figures  showed  that  the  spectra  of  normalized  power  for  perturbed  data  with  density 
fluctuations  displayed  finer  scale  structures.  In  the  volume  plot,  a  distinguished  drop-off 
of  normalized  power  beyond  a  wavenumber  of  approximately  0.1  m"'  was  observed  for 
the  perturbed  data  with  density  fluctuations. 

Figures  27  and  28  presented  the  results  for  the  broadband  signal  for  the  interface 
and  volume  reverberation.  Large-scale  structures  were  observed  in  both.  The  perturbed 
data  with  density  fluctuations  did  not  significantly  change  the  results  from  the  perturbed 
data  without  density  fluctuations. 
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Source  depth:  48m,  Receiver  depth:  40m,  Frequency:  250Hz 


Figure  25.  Normalized  Power  Spectrum  of  Interface  Reverberation  (CW) 


Wavenumber,  (m"') 


Figure  26.  Normalized  Power  Spectrum  of  Volume  Reverberation  (CW) 
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Source  depth:  48m,  Receiver  depth:  40m 


Figure  27.  Normalized  Power  Spectrum  of  Interface  Reverberation  (Broadband) 


Source  depth:  48m,  Receiver  depth:  40m 


I  —  Perturbed  Power  without  Density  Fluctuations  I 

I  —  Perturbed  Power  with  Density  Fluctuations  | 


Wavenumber,  (m"') 


Figure  28.  Normalized  Power  Spectrum  of  Volume  Reverberation  (Broadband) 
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2.  The  Analysis  of  Power  Ratio  Spectral  Density  (PRSD) 

In  order  to  further  examine  speetral  content,  the  power  ratio  spectral  density  was 
implemented.  It  was  via  the  use  of  discrete  Fourier  transforms  (DFT)  of  the  magnitude 
square  of  the  ratio  of  the  reverberation.  It  was  then  defined  as'^'^ 


PRSD  =  DFT  { 


P  perturbed 

2 

P  unperturbed 

f" 

(3.11) 


For  the  CW  signal  and  perturbed  data  with  and  without  density  fluctuations, 
Figures  29  and  30  portrayed  the  normalized  spectral  components  of  power  ratio  for  the 
interface  and  volume  reverberation,  respectively.  With  density  fluctuations  in  the 
perturbed  data,  the  level  of  normalized  power  appeared  to  be  less  than  the  data  without 
density  fluctuations.  Additionally,  at  greater  wavenumbers,  the  volume  reverberation  data 
displayed  a  downward  structure  while  the  interface  showed  a  fiat-like  structure. 

Lastly,  the  computation  of  power  spectral  density  for  the  broadband  signal  was 
performed  based  on  the  time-series  data.  The  results  for  the  interface  and  volume  were 
shown  in  Figures  31  and  32.  The  observation  for  these  plots  was  agreed  with  the  CW 
analysis  for  the  ratio  between  the  perturbed  and  unperturbed  data  sets. 
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(|Pertijrbed|^/|Unpertijtbed|^)  Source  depth:  48m,  Receiver  depth:  40m,  Frequency:  250Hz 


I  Perturbed  Data  without  Density  Fluctuations  I 
I  —  Perturbed  Data  with  Density  Fluctuations  | 


10^  10"^  10  '  10° 

Wavenumber,  (m”’) 


Figure  29.  Normalized  Spectrum  of  Interface  Power  Ratio  (CW) 


(|Per1urbed|^/|Unper1urbed|^)  Source  depth:  48m,  Receiver  depth:  40m,  Frequency:  250Hz 


10^  10”^  10  '  10® 

Wavenumber,  (m”') 


Figure  30.  Normalized  Spectrum  of  Volume  Power  Ratio  (CW) 


43 


(|Per1urbed|^/|Unper1urbecl|^)  Source  depth:  48m,  Receiver  depth:  40m 


I  —  Perturbed  Data  without  Density  Fluctuations  I 
I  —  Perturbed  Data  with  Density  Fluctuations  | 


10^  10"^  10  '  10° 

Wavenumber,  (m"') 


Figure  3 1 .  Normalized  Spectrum  of  Interface  Power  Ratio  (Broadband) 


(|Petturbed|^/|Unperturbed|^)  Source  depth:  48m,  Receiver  depth:  40m 


10’^  lO”^  10’’  10® 

Wavenumber,  (m"’) 


Figure  32.  Normalized  Spectrum  of  Volume  Power  Ratio  (Broadband) 
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IV.  SUMMARY 


The  focus  of  this  thesis  has  been  the  use  of  the  improved  MMPE  model  to 
investigate  the  influence  of  density  fluctuations  for  both  interface  and  volume 
reverberation.  Experiments  were  conducted  using  the  same  reverberation  geometry  and 
environmental  parameters  as  defined  in  the  previous  work.  A  vertical  line  array  (VEA) 
with  16-elment  was  chosen  to  provide  the  reverberation  measurements  and  computations. 
The  source  was  located  at  depth  of  48  m  with  all  16  elements  receiving  the  reverberation. 
In  addition,  a  relative  mean  square  (rms)  value  of  1  m  interface  roughness  and  a  15  m/s 
volume  sound  speed  fluctuations  were  used.  To  compare  the  perturbed  data  with  and 
without  density  fluctuations  in  the  environment,  several  analyses  (reverberation  pressure 
levels,  vertical  correlation,  peak  correlation  and  spectral  characteristics)  were  performed 
for  both  CW  and  broadband  signals. 

The  first  method  of  analysis  was  to  observe  the  structure  of  mean  reverberation 
pressure  levels  (RPE).  Plots  of  reverberation  loss  versus  range  and  time  were 
accomplished  for  CW  and  broadband  analysis,  respectively.  We  noted  that  the  structure 
of  mean  reverberation  pressure  levels  remained  statistically  uniform  over  region  of  study 
assuming  the  dominant  scattering  mechanism  was  small-scale  Bragg  scatter. 
Additionally,  the  influence  of  bottom  volume  density  perturbations  was  to  reduce  later 
levels  relative  to  earlier  levels  but  did  not  appreciably  affect  structure.  This  was  probably 
a  result  of  direct  correlation  between  volume  sound  speed  and  density  fluctuations. 
Noticed  that  the  CW  analysis  was  unable  to  capture  coherent  structure  of  volume  RPE 
due  to  the  inability  to  resolve  multipath  influence. 

Next  we  considered  the  vertical  correlation  of  interface/volume  reverberation  loss. 
Eor  CW  analysis,  the  structures  of  perturbed  data  with  and  without  density  fluctuations 
were  similar  for  both  cases  of  the  interface  and  volume  reverberation.  Eor  the  interface 
reverberation  range  correlation,  the  perturbed  data  with  density  fluctuations  revealed 
stronger  correlation  throughout  depth  for  the  range  lag  between  0.3  and  0.5  km.  This 
stronger  correlation  could  also  be  found  in  the  perturbed  data  with  density  fluctuations 
for  the  range  lag  between  2.7  and  3  km  in  the  analysis  of  volume  reverberation. 
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Comparison  of  perturbed  data  with  and  without  density  fluctuations  for  broadband 
analysis  showed  that  the  difference  between  them  was  very  small.  This  tiny  difference 
applied  to  both  the  interface  and  volume  reverberation  temporal  correlation.  Furthermore, 
the  finer  scale  structures  were  also  noted  in  both  the  interface  and  volume  reverberation 
temporal  correlation. 

Then  in  order  to  observe  the  rate  of  change  for  the  signal  decorrelated  over  depth, 
the  results  on  the  perturbed  data  computed  with  and  without  density  fluctuations  of  peak 
vertical  correlation  were  extracted  from  the  previous  analysis  of  vertical  correlation  for 
the  interface  and  volume  relative  to  a  source  depth  of  48  m.  Because  of  the  inability  to 
resolve  multipath  influence,  peak  vertical  correlation  analysis  only  valid  for  broadband 
pulse  calculations.  According  to  the  results,  the  peak  vertical  correlation  analysis 
suggested  that  volume  reverberation  may  decorrelate  across  vertical  array  more  rapidly 
than  interface  reverberation.  This  was  presumably  due  to  multi-point/multi-depth  scatter 
contributions  of  the  volume  producing  more  vertical  structure  than  the  interface. 

Lastly,  to  extract  the  spectral  components  of  the  reverberation  data  for  both  the 
CW  and  broadband  signals,  the  analysis  was  performed  by  using  the  methods  of  signal 
processing,  such  as  the  power  spectral  density  and  power  ratio  spectral  density  .  Spectral 
analysis  of  both  CW  and  broadband  pulse  calculations  suggested  that  the  response  of 
interface  reverberation  was  flatter  with  a  slope  on  the  order  of -0.125  for  both  CW  and 
broadband  data.  Volume  response  displayed  a  -0.75  slope  for  CW  and  around  -0.25 
slope  for  broadband.  There  was  no  clear  connection  between  spectral  characteristics  of 
predicted  signal  and  environmental  perturbations. 

With  the  conclusion  of  this  thesis,  recommendations  for  future  work  are  as 
follows: 

•  Perform  short-range  statistical  analysis  to  narrow  down  relationship  between 
signal  structure  and  environmental  structure.  Attempt  to  interpret  influence  of 
multipath  effects  on  long-range  structures. 

•  Incorporate  environmental  measurements  from  ASIAEX  data  and  begin 
prediction  analysis. 
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•  Begin  data  processing  of  ASIAEX  data  and  perform  data/model  comparisons. 

•  Incorporate  rough  surface  scatter  into  propagation  model  and  investigate 
influence  on  various  RPL  predictions. 
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